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ABSTRACT 

We present the results of our optical monitoring of the BL Lac object 
S5 0716+714 on seven nights in 2006 December. The monitoring was carried 
out simultaneously at three optical wavelengths with a novel photometric sys- 
tem. The object did not show large-amplitude internight variations during this 
period. Intranight variations were observed on four nights and probably on one 
more. Strong bluer-when-brighter chromatism was detected on both intranight 
and internight timescales. The intranight variation amplitude decreases in the 
wavelength sequence of B\ R', and V. Cross correlation analyses revealed that 
the variability at the B' and V bands lead that at the R' band by about 30 
minutes on one night. 
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Subject headings: BL Lacertae objects: individual (S5 0716+714) — galaxies: 
active — galaxies: photometry 



INTRODUCTION 



Blazars represent the most violently variable objects among all active galactic nuclei 
(AGNs). They show rapid and strong variability, high and variable polarization, and a 
non-thermal continuum. The variable continuum is believed to come from near the base of 
the relativistic jet closely aligned with our light of sight. The jet is probably powered and 
accelerated by a rotating supermassive black hole surrounded by an accretion disk. A blazar 
can either be termed as a flat-spectrum radio quasar (FSRQ) or a BL Lac object, depending 
on whether or not it shows strong emission lines in its spectrum. 

The spectral energy distribution (SED) of blazars has a two-bumped structure: The low- 
frequency bump ranges from radio to UV or X-ray frequencies and is believed to come from 
the synchrotron emission from the relativistic electrons in the jet. The high-frequency bump, 
from X-rays to 7-rays, is usually interpreted as emission from inverse Compton scattering 
of the low frequency radiation by the same ensemb le of relativistic electrons responsible for 
the synchrotron emission (for a recent review, see iBottcherl 120071 ). According to the peak 
frequency of their synchrotron emission, BL Lac objects can be further divided into low- 
frequency-peaked BL Lacs (LBL) and high-frequency-peaked BL Lacs (HBL). The former 
subclass has the synchrotron peak at IR-optical wavelengths while the latter subclass has 
the peak in the UV or soft X-ray band. 



The BL Lac object S5 0716+714, probably at a redshift of 0.31+0.08 (INilsson et al. 



2008 ). is one of the best studie d blazars. It is highly variab 
1996 ) and 7-ray wavelengths (jLin et al. 1995 : Chen et aL 



tected very high energy 7-ray emission from this object (I Ander hub et al.ll2009l). Its SED has 



e from radio to X-ray (jWagner et al. 



20081). Most recentl v. MAGIC de- 



been obseryed and studied by a number of multi-wavel e ngth campaigns (e.g 



199 



20081: 



Tadia 



'erri et al. 2003: Ferrero et al. 2006; Foschini et al. 2006 



Giommi et al. 



Ostorero et al. 



20061 : Ivillata et al.lboosi ). The majority of the observations revealed a two-bumped structure 
of the SED with a concave shape at 2 - 3 keV. Furthermore, it was found that its emission 
from radio to soft X-rays, likely due to the synchrotron process, could be highly variable, 
while the hard X-ray (from 3 to 10 keV) flux, most probably dominated by inverse Compton 
emission, usually remained constant]^ A model with two synchrotron self-Compton (SSC) 
components was proposed to explain its SED changes, with one component highly variable 



^An exception was found by Ferrero et al. ( 2006I ). who detected its delayed hard X-ray variation as well. 
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and the other constant (jChen et aLll2008l : iGiommi et aLl 120081 ) 



In the optical regime, S5 0716+714 is one of the brightest blaz ars. It has a duty cycle 



close to 1, which means that it almost n ever stops its variation (e.g.. IWagner fc Witzellll995 



timescales by a number of authors (IGhisellini et al. 



2000 



Nesci et al 



2002 



2005 



Stalin et al.ll2009t IChandra et al.l 120111). Its optical variabil i ty has been studied on different 



Qian et al. 



2002 



1997 



.-laiteri et al. 



Sagar et al. 



2003 



1999 



Xie et al. 



rvii 



2004 



ata et al. 



Wu et al. 



2005 



2007al: 'Monta gni et al. 



2006; Stalin et al. 



2006 



Pollock et al. 



Sasada et al.l 



Carini et al 



2008j : iZhang et al.ll2008i: iStahn et al.ll2009[ IPoon et al.lbood i lzhang et al 



20111 



Chandra et al 



maghr ^ (ISagar et al 



1999 



200 



. Gupta et a 


2008 


Zhang et al. 


2010a 



2011). Its variation rate can be faster than 0.10-0.12 



Villata et al. 2000; Wu et al 



2005 



Montagni et al.l 120061 ). A 



maximum rate of 0.16 maghr"^ was reported by iNesci et al.l (120021). The variation timescal e 
can be as shor t as 15 minutes ( Sasada et al. 20081 ; Rani et al. 2010a; Chandra et al. 2011 ). 
Most recently, iFan et al.l ( 1201 ll ) even reported a 0.611 mag variation over 3.6 minutes in this 
object. 

The color or spectral behavior of blazars is a subject o f much debates. So me au- 



thors found a bluer- when-brighter (BWB) chromatism (e.g., IVagnetti et al.l l2003l ). some 
others claimed the opposite behavior, namely, a redder- when-brigh t er (R WB) trend (e.g.. 



Ramirez et al.ll2004( ). or no clear tendency ( e.g.. iBottcher et al 



BL Lacs are BWB while FSRQs are RWB dPan fc LinI 



2000 



2007 



20091). It appears that 



Gu et al. 2006; Hu et al. 2006 



Rani et al.lboiObhl^ A certain object rn ay display different trends in dif ferent variation modes 



Wu et al. 



2005; Poon et al 



n— i 20091 ) or on different timescales (e.g.. IGhisellini et al.lll997l ; 

Raiteri et al. 20031 ) . One reason for these divergences is that almost all previous observations 
were done quasi-simultaneously at multiple wavelengths. When the exposure was switched 
from one wavelength to the other, the brightness of the source might have changed. So the 
quasi-simultaneous observations cannot obtain real color, and the color behavior drawn from 
these observations should be taken with caution. 

We monitored S 5 0716-1-714 for seven nights in 2006 December with a novel optical 
photometric system ( Wu et al. 2007al lbl). This system enables simultaneous photometry at 
multiple optical wavelengths, and thus has the advantages to accurately trace the color 
change of blazars and to increase the temporal resolution at individual wavelengths. The 
resulting data were used to study the properties of the short-term variability of S5 0716+714, 
to analyze its color behavior during variations, and to search for the correlation and time lags 
between the variations at various optical wavelengths. Section 2 describes our observations 
and data reduction procedures. Section 3 presents the light curves and the results of the 



^However, 



Gu fc Ai 



(|201l[ ) found that only one FSRQ is RWB in a sample of 29 SDSS FSRQs. 
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color and correlation analyses. The conclusions and discussions are given in Section 4. 



OBSERVATIONS AND DATA REDUCTION 



The monitoring of S5 0716+714 was carried out at Xinglong Station of the National 
Astronomical Observatories of China in 2006 December with a novel photometric system, 
which consists of a 60/90 cm Schmidt telescope, an objective prism, a CCD camera, and 
a multi-peak interference filter (MP IF). Beams of light in multiple passbands are refracted 
differentially by the objective prism, reflected by the spherical main mirror, pass through 
the MPIF simultaneously, and are focused on the CCD separately. So we have multiple 
"images" on the frame for each object. The left three passbands of the MPIF are close to 
the broadband B, V, and R, so we designate them as B', V, and R'. A fourth passband 
beyond 9400 A is neglected, b ecause the CCD h as a very low response at those wavelengths. 
For details of the system, see IWu et al.l ( l2007al Jbl). 



A change to the system is that we now use a new 4096 x 4096 E2V CCD. Compared to 
the old one, it has a smaller pixel size of 12 /jm and a higher spatial resolution of 1.3 "pixel"^. 
When observing blazars, we read out only the central 512 x 512 pixels as a frame, which has 
a field of view of about 11' x 11'. The key advantage of the new CCD is its high quantum 
efficiency at blue wavelengths (92.2% at 4000 A). So we now have usef ul data in one ni ore 
band, the B' band, in addition to the V and R' bands for the old system ( Wu et al. 2007al lbl). 



The monitoring covered the period from 2006 December 18 to 27 (from JD 2,454,088 
to 2,454,097), with no observation on December 19, 25, and 26 due to bad weather. The 
exposure times are all 150 seconds. Taking into account the readout time of about 6 seconds, 
we achieved a temporal resolution of less than three minutes in three passbands. The average 
FWHMs of the stellar images are 3'.'1, S'.'l, 3'.'6, 4'.'3, 3'.'4, 3'.'7, and 5'.'4 on JDs 2,454,088, 
2,454,090, 2,454,091, 2,454,092, 2,454,093, 2,454,094, and 2,454,097, respectively. 

An example frame is shown in Figure [H We rotated the objective prism (and thus the 
dispersion direction) by a small angle from the standard vertical position. As the result, 
some overlaps of images were resolved, especially the potential overlap of the B' image of 
S5 0716+714 with the i?' ima ge of a star below (for a comparison, see Fig. [T]in this paper 
and Fig. 2 of IWu et al.ll2007ah . 



The data reduction included bias subtraction, flat-flelding, extraction of instrumental 
aperture r nagnitu de, and flux calibratior i. We defined a specially shaped aperture, as in 
Wu et al.l (l2007a] ). but with an update. In lWu et al.l (l2007al ). the traditional circular aperture 
was equally split into two semicircles and a rectangle was inserted between them. Then an 
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elongated aperture was obtained for the elongated images. For the inclined and elongated 
stellar images in this paper, we replaced the rectangle with a parallelogram. The chords of the 
two semicircles were attached to the top and bottom sides of the parallelogram, respectively, 
while the left and right sides of the parallelogram were set parallel to the dispersion line. 
Different parallelogram heights were adopted for the apertures in different passbands. They 
are 11, 7, and 4 pixels for the B', V, and R' images, respectively. The corresponding 
parallelogram widths (or the diameters of the semicircles) were set to be 5, 6, and 6 pixels, 
respectively. All sky annuli were set to have inner and outer radii of 11 and 15 pixels, but an 
area within 6 pixels in right ascension to the central axis of the aperture was excluded in order 
to e liminate the possible contamination of the background by the apertures. Stars 1,2, and 
4 in IVillata et al.l (Il998[ ) were used as reference stars. Their average brightness was utilized 
to calibrate the flux of S5 0716+714. Star 6 is slightly brighter than S5 0716 +714 and was 



used as a check star (for a reasonable selection of reference and check stars, see iHowell et al. 



19881 ). The above aperture sizes in the three bands were chosen based on several trials to 



minimize the scatter s in the d i fferenti al magnitudes between the check star and the three 
reference stars. As in lWu et al.l ( l2007al ). the standard B', V, and R' magnitudes of these four 
stars were mimicked by their broadband B, V, and R magnitudes. There should be some 
differences between these two magnitude systems, but they are tiny ( IWu et al.ll2007al ). and 
are irrelevant for the results of our variability study presented here. The results are presented 
in Table 1. The columns are observation date and time (UT), Julian Date, exposure time, 
and magnitudes and errors of S5 0716+714 and differential magnitudes of star 6 in the B', 
V, and R' bands, respectively. 



RESULTS 



3.1. Light Curves 



Figure [2] shows the total light curves of S5 0716+714 in the B' , V', and R' bands. The 
B' and R' band light curves are shifted by 0.2 and —0.5 magnitudes, respectively, in order 
to make the variations clearly visible. The nightly-averaged brightness of this object was 
quite stable in this period and varies by less than 0.2 mags. On the other hand, intranight 
variations can still be observed on most nights by visual inspection. 

We then performed a quantitative assessment on whether there are variatio ns on these 
seven night s. Here we used t he one-way analysis of variance (ANOVA) proposed by lde Diego et al, 
(Il998[ ) and Ide Diego (120101 ). It is proved to be a powerful and robust estimator for varia- 
tions ( de Diego[2010 ). In brief, the observations on each night were divided into groups of 
five consecutive data points. The variance of the means of each group about the overall 
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mean was computed, as was the mean of the variances or dispersions (residuals) within each 
group. Then the ratio between these two variances was calculated and multiplied by five, 
the number of observations in each group. The value obtained behaves as the F statistic. 
For a certain significance level, if F exceeds the critical value, the null hypothesis that there 
is no variation will be rejected. The results are presented in Table 2. Column 1 is the Julian 
Date. Columns 2 and 3 are the degrees of freedom for the groups and residuals, respectively. 
Columns 4, 5, and 6 are F values in the B', V, and R' bands, respectively. Column 7 is the 
critical value at the 99% significance level at the specified two degrees of freedom of Columns 
2 and 3. Column 8 indicates whether there is variation. As can be seen, there are variations 
on JDs 2,454,088, 2,454,090, 2,454,092, and 2,454,097. The F^', Fy, and Fr/ are all greater 
than the critical values on the four nights. On one more night, JD 2,454,091, there may be 
variation, too. The Fyi and Fr/ are greater than the critical value, whereas the F^/ is slightly 
smaller than it on that night. On the two remaining nights, JDs 2,454,093 and 2,454,094, 
there is no variation. 

Figures [3l HJ |5l and |6] show the intranight light curves of S5 0716+714 on the four 
nights with variations, JDs 2,454,088, 2,454,090, 2,454,092, and 2,454,097, respectively. The 
associated small panels give the differential light curves of the check star (the means were 
set to 0). For t hese four nights, we ca lculated the intranight variation amplitudes with 
the definition of iHeidt fc Wagnerl ( 19961 ). A = ^J (mmax — ^rain)'^ — 2o"^, where mmax and 
mmin are the faintest and brightest magnitudes during any given night, respectively and 
a is the standard deviation of the 'variation' of the check star. The last point in the R'- 
band light curve on JD 2,454,088 and the point at (0.26764, 14.034) in the S'-band light 
curve on JD 2,454,097 are likely to be spurious measurements (see Figs. |3]and[6]) and were 
discarded before the calculations. The results are listed in Table 3. It can be seen that 
the i?'-band variations have the largest amplitudes and the V- and i?'-band amplitudes 
are similar, with the latter lightly larger than the former on three of the four nights. In 
the shock-in-jet model, the variation am plitude decreases towards longer wavelengths . This 



has been observed in S5 07 16+714 (e.g.. lWu et al. 



and in other blazars (e.g., iNesci et al.l Il998t iFan fc LinI l2000t iDai et al 



2007a 



Poon et al 



2009; 



Fan et al.l 120091) 



201 If ). However, 



the amplitudes may also be related to the relative position of the frequency at which the 
observations were made and the peak frequency of the synchrotron compone nt in the SED 



of blazars. The latter frequency itself may also shift from time to time (e.g., iGiommi et al. 



19991 . l2000t iTavecchio et al.l 1200 ll : IZhang et al.ll2010bl ). Taking these facts into account, our 
results that the i?'-band ampli tudes are s lightly larg er than the y^-b and ones may not be 
unreasonable. As mentioned in Iwu et al.l (2007a) and Dai et al. (2011), the larger variation 
amplitude at higher frequency can partly explain the BWB chromatism observed in blazars. 
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3.2. Color-Magnitude Correlation 



Controversy exists in tlie studies of tlie spectral or color behavior of blazars, as men- 
tioned in §1. For our data on the four nights with variations, JDs 2,454,088, 2,454,090, 
2,454,092, and 2,454,097, the color indices B' - V and B' - R' were calculated and plotted 
against the B' magnitudes in the left and right panels of Figure [TJ respectively. Different 
symbols mark data on different nights. They were fitted linearly and separately for each 
singl e night. We used the BCE S regression (bivariate correlated errors and intrinsic scat- 



ters, lAkritas fc Bershadyl Il996l ) to do the fitting. The BCES regression has the advantage 
of taking into consideration not only the measurement errors, but also the intrinsic scatters. 
The regression lines are plotted in Figure [7] with different styles for different nights. The 
thick solid lines are the fittings to all points. The correlation coefficients are given at the 
upper-left corners. It can be seen that the object shows strong BWB chromatisms on both 
intran ight and intern i ght timescales. This result is in agreement with several previous re- 



sults flWu et al.1 12005 



2007a: Hao et al.l 2010l: IChandra et al. 2011) and with the short-term 



behavior obtained by iGhisellini et al.l (Il997( ) and IZhang et al.l (|2008[ ). The much larger cor- 
relation coefficient in the right panel indicates a stronger correlation between B' — E! and 



B' than between B' — V and B' . A similar result has been reported by lRaiteri et al.l ( l2003l ). 
who found a stronger BWB trends with the B — I color than with the B — R color. 



In addition to the flatter-when-brighter trend during rapid flares, IGhisellini et al.l (Il997| ) 
have reported that the correlation was rather insensitive for the long-terra variations of 
S5 0716+714. A similar phenomenon was also recognized by iRaiteri et al.l (|2003[ ) in the 
same object. Our monitoring lasted for only ten days, so we cannot address the long-term 
color behavior of this object. 



3.3. Color Evolution During Flares 



In a variety of blazar emission models, one expects an interplay between gradual particle 
acceleration, radiative cooling and escape. This will lead to a loop-like path of the blazar's 
state in a color-magnitude (or spectral index-flux) diagram. The direction of this spectral 
hysteresis can be either clockwise or anticlo ckwise, depending on the relative values of the 
acceleration, cooling and escape timescales (IChiaberge fc Ghisellinilll999l : lDermerlll998l ). as 
well as the frequency at which the observation is made relative to the peak frequency of the 
synchrotron component in the SED of the blazar (see F igs. 3 and 4 in 



Kirk et al. 1998). The 



loop path has been i 


requently reported in X- 


1996; 


Kataoka et al. 


1999; 


Zhane: et al. 


1999[ 



and, in one case, in the infrared ( IGear et al.l Il986l ). In the optical regime. 



1993 



Takahashi et a. 



Xilouris et al. 
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( 120061 ) have reported a similar pattern, and IWu et al.l ( l2007a( ) showed a clockwise loop with 



the nightly means of magnitudes and colors. Such hysteresis behavior will inevitably translate 
into inter-band time lags. 

In our monitoring, there are some mini-flares. The one with the largest amplitude 
occurred from 0.245 to 0.292 on JDs 2,454,088 (see Fig. [3]). We checked the color evolution 
during that flare. Figure |8] displays the result. The numbers associated with the points 
signify the evolutionary sequence, with the beginning and end points marked by the filled 
circle and square, respectively. There is a hint of a clockwise loop when we check only 
the points, but the relatively large errors prevent us from drawing a convincing conclusion. 
Therefore, in order to observe an optical hysteresis loop on the color-magnitude diagram, 
more accurate observations are needed, preferentially during a fiare with a larger amplitude 
than that of the current one. High temporal resolution may also be crucial. 



3.4. Cross Correlation Analyses and Time Lags 

We then performed the inter-band correlation analyses and searched for the possible 
inter-band time delay. These are carried out on the intranight light curves for the four 
nights with variations, JDs 2,454,088, 2,454,090, 2,454,092 and 2,454,097. 



We at first used the z-transformed discrete correlation functions (ZDCFs; [Alexander 



19971 ) to search for the B'-V, V'-R' and B' R' correlations . The ZDCF differs from the 



discrete correlation function (DCF) of lEdelson fc KrolikI ( 1l988l ) in that it bins the data points 



into equal population bins and uses Fisher's z-transform to stabilize the highly skewed dis- 
tribution of the correlation coefficient. It is much more efficient than the DCF in uncov- 
ering correlations involving the variability timescale, and deals with under-sampled light 
curves better than both the DCF and the interpolated cross-correlation function (ICCF; 



Gaskell fc PetersonI 119871 ). For each correlation, a Gaussian fitting (GF) was made to the 
central ZDCF points with ZDCF value greater than ~ 75% of the peak value. The time 
where the Gaussian profile peaks denotes the lag between the correlated light curves. 

The results are listed in Table 4 and displayed in Figures |9| [T0 | [TT | and[T2]for JDs 2,454,088, 
2,454,090, 2,454,092, and 2,454,097, respectively. In each panel, Julian Date and the corre- 
lated passbands {B'-V , V'-R', or B'-R') are given at the upper-left corner. The peak of the 
Gaussian profile (the dashed line) is marked with a vertical dotted line. The corresponding 
lag is given at the upper-right corner. A positive lag means that the former variation leads 
the latter one. 



One problem in the GF is that it usually gives an unreliable and under-estimated error 



- 9 - 



for the lag due at least to the statistically non-independent ZDCF points. Then, although 
the lags obtained from the ZDCF+GF method generally have significance greater than 3cr, 
they should be taken with caution. 



Another way to measure the lags and their errors is the ICCF (IGaskell fc Petersonlll987l ). 
The error was estimated with a model-independent Monte Carlo method, and the lag was 
taken as the centroid of the cross-correlation functions that were obtained with a large num- 
ber of independent Monte Carlo realizations. Thi s is the flux-ran d omiza tion / random- subset 
selecti on (FR/RSS) approach, as described by iPeterson et al.l (Il998l ) and iPeterson et al. 
( 120041 ). In the current work, five thousand independent Monte Carlo realizations were per- 
formed on the light curves on JDs 2,454,088, 2,454,090, 2,454,092, and 2,454,097, and the 
lags and errors were obtained and are presented in Table 4. One can see from the table 
that although the ZDCF+GF and FR/RSS lags are usually consistent with each other, the 
FR/RSS errors are much larger than the corresponding ZDCF-I-GF errors. The FR/RSS 
lags have significance lower than 3a except for the V — R' and B' — R' lags on JD 2,454,090. 
On that night, the variability at the B' and V' bands lead that at the R' band by about 30 
minutes. 



A few optical lags have been reported in several blazars before. For example, iRomero et al 

( 120001 ) derived a ~ 4 minute lag of the R band variations relative to the V band variations for 



PKS 0537-441 (see their Fig. 5). BL Lac was observed t o have its / band varia tions delayed 
by ~ 0.2 hours with respect to its B band variations ( iPapadakis et al.ll2003l ). The same 
authors also reported a lag of 0.9 ±1.0 hours of the / band variations rela tive to the B band 
variations for S4 0954+658 JPapadakis et al.l 120041 ). For S5 0716+714, loian et al.l J2OO0I ) 
found an upper limit o f 6 mi nutes for the time lag between the V and / band variations. 
Similarly, IVillata et al.l (j2000f ) reported an upper li mit of 10 mi n utes o n the delay of the / 
band variations relative to the B band variations. IStalin et al.l ( l2006l ) also presented time 
lags of 6 and 13 minutes for the V and R variations on two nigh ts. But they also c autioned 
about the low significance level of their results. Most recently, iPoon et al.l (l2009l) f ound a 
possible lag of about 11 minutes between B and / band variations, and IZhang et al.l ( l2010al ) 
reported possible lags of a few minutes between low and high optical frequencies. These 
previous lags are either upper limits or cautioned by the authors to have low confidence lev- 
els. In the current work, however, the variability at short wavelengths were observed to lead 
that at long wavelengths on at least one night. This positive detection of inter-band time 
delays presumably has taken advantage of the simultaneous measurements at three optical 
wavelengths of our novel photometric system. 
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4. CONCLUSIONS AND DISCUSSIONS 
4.1. Conclusions 

We monitored S5 0716+714 for 7 nights with our novel photometric system in 2006 De- 
cember. The photometry was carried out simultaneously at three optical wavelengths. The 
object did not show large-amplitude internight variations during this period, but displayed 
intranight variations on four nights and probable variations on one more. Strong BWB chro- 
matism was detected on both intranight and internight timescales. For the four nights with 
variations, the intranight variation amplitude decreases in the wavelength sequence of B', 
R', and V. Cross correlation analyses revealed that the variations at the B' and V bands 
led that at the R' band by about 30 minutes on one night. 



4.2. Discussion on Color Behavior 



S5 0716+714 showed a strong BWB chromatism during this period, either when the data 
are considered together or separately for individual nights. Despite the strong correlation, 
the points are in fact quite diffuse on the color- magnitude map (Fig. [7j). Data on different 
nights occupy quite different regions, and follo w distinct color- magnitude correlations. This 
disp ersion is a l so ev ident in IWu et al.l (j2007al ) (their Fig. 9). As mentioned in that paper 
and iDai et al.l (120111 ). differences in variation amplitudes and paces (the lag) can both lead 
to the observed BWB chromatism. The variation amplitudes and the possible lags presented 
in this paper vary from night to night. So the dispersion on the color- magnitude map may 
be linked to the dispersion on the intranight variation amplitudes and/or the possible lags. 

For the spectral or color behavior of blazars, different authors have drawn different 
conclusions. However, their results on the color-magnitude correlation of blazars may be 
biased by many aspects. Here we discuss some of them. 



1. Non-simultaneous measurements at different wavelengths. Strictly speaking, the spec- 
tral or color index should be calculated with simultaneously measured fluxes or mag- 
nitudes at different wavelengths. However, most observations can only obtain quasi- 
simultaneous measurements. For blazars which do not show significant variations on 
very short timescales (< 0.5 hour), this may not be a problem. However, for those 
showing very fast variations, like S5 0716+714, as mentioned in the first section, this 
may be a serious problem. The intrinsic brightness of the object may change measur- 
ably when two observations are made successively at separated wavelengths. These two 
observations will result in a wrong spectral or color index. Although an interpolation 
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can ameliorate this problem to some extent, the irregularity of the variability of blazars 
makes the interpolation usually unpractical. Our novel photometric system with the 
MPIF has the ability to take measurements simultaneously at three wavelengths. It 
can thus trace color changes accurately. 

2. Observations taken at multiple telescopes and/or sites. There are likely systematic 
differences between data o btained with different telescopes. The difference can be 0.05 
mag or even higher (e.g., iGhisellini et al.l 119971 1. Then the color index will have an 
error close to 0.07-0.1 mag, a value comparable to the scale of the color change. This 
error will definitely weaken the color-magnitude correlation. 

3. Inclusion of other components of variation. In addition to the shock-in-jet component, 
other mechanisms may be involved and dominant for an episode in blazar va riability, 
such as the ge ometric effects suggested for S4 0954-1-658 (IWagner et al.l 119931 ) and for 
S5 0716-1-714 ( IWu et al.ll2005l ). These effects lead to achromatic variations. In another 
BL Lac object, OJ 28 7, a huge achromatic outburst during 1994-96 was reported by 



Sillanpaa et al.l (Il996l ). However, this outburst was dominated by a component pre- 



sumably resulting from the interaction of the primary accretion d i sk and the secondary 



Liu & Wu 


2002 




Valtonen et al. 


20081) 


showed a BWB trend ( 


Wu et al. 


2006[ 



20081) . When the blazar variability was isolated, it 



4. Utilization of the averaged fluxes. If, for a blazar, the difference between variation 
amplitudes at different wavelengths is very small, and the spectral change is mainly 
resulting from time lag, we will not derive a spectral change when averaging the fluxes 
within a timescale far l onger than the tim e lag. Therefore, it is not surprising that 



Sagar et al.l (119991 ) and IStalin et al.l (120061 ) claimed no spectral change, because they 



used the averaged flux or magnitude to derive the spectral or color index. However, 
when the amplitude difference in a blazar dominates the color change, the color change 
may also be detected by using averaged fluxes. 

5. Utilization of closely separated frequencies. The B — V and V — R color indices change 
m uch less than the B — R, V — I, and B — I color indices do. This has been illustrated 
by iRaiteri et al.l (120031 ) and this work. 



There may be other barriers for the detection of the color or spectral changes and the 
color-magnitude correlations. These are the five major issues. Some of them may combine to 
function in the observations. For exampl e, a campaign involviri g sev eral telescope s usua lly 
suffer from the first two, as in the case of IGhisellini et al.l (119971 ) and iRaiteri et al.l (120031 ). 
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Our photometric system has the advantage of simultaneous measurements at three opti- 
cal wavelengths, and suffers none of the above 5 barriers in the detection of the color change. 
So the resulting color behavior should have higher confidence level than those from previous 
non- or quasi- simultaneous observations. 



4.3. Discussion on Time Lag 



A few intrinsic or extrinsic models have been proposed to explain the variability of 
blazars. The most popular is the shock-in-jet or internal shock model, in which shocks 
propagate down the relativistic jet, accelerating particles and/or compressing magnetic fields, 
leading to the observed variability. Depending on the balance between escape, acceleration, 
and cooling of the electrons wit h different energy, either soft (low energy) or hard (high 
energy) lags are expected (e.g., iKirk et al.l Il998t ISokolov et al.l l2004l : iMarscherl |2006| ). In 
fact, time lags between variations at much separated wavelengths are well known. For 
example, the soft X-ray variatio ns were observed to lag the hard X-ray s by about 1 hour i n 
Mrk 421 jTakahashi et al.lll996h and by about 1000 s in S5 0716+714 Jzhang et al.lboiObI ). 
In PKS 2155 -304, the X-ray fl are appears to lead the EUV and UV fluxes by 1 and 2 days, 
respectively ( lUrry et al.l 119971 ). Sometimes, the lags may show as the offsets between the 
X-ray knots and the radio or optical one s in the relativis tic jets of radio galaxies, with the 
X-ray knots being closer to the core (e.g. 



Bai fc Led 120031 ) 



I n optical regime, se veral authors claimed no detection of inter-band delays (e.g. lHao et al 



2010l : ICarini et al.l 1201 ll ). Some other claimed the opposite, as mentioned in §3.4. Four key 
parameters probably determine whether the optical time lags can be detected: (1) wave- 
length separation, (2) variation amplitude, (3) temporal resolution, and (4) measurement 
accuracy. For the first factor, the larger the wavelength separation, the easier the lag will 
be detected. Then a.B — RotB — I correlation should be performed rather than a. B — V 
or V — R correlation. For the second, blazars with high duty cycles and large variation 
amplitudes should be observed for the search of the lag. For the third, a temporal resolution 
of less than 5 minutes is favorable, because most reported time lags are less than 10 minutes, 
as mentioned above. For the fourth, bright sources should be monitored. S5 0716+714 is one 
of the brightest blazars in the optical regime. It shows a high variation rate and has a high 
duty cycle, as mentioned in §1. This makes it probably the best candidate to be monitored 
for the search of the optical lags. 

Our photometric system has the advantage of simultaneous measurements at multiple 
wavelengths and increases the temporal resolution significantly in each single wavelengths. 
With this system, we detected the time delays on one night. On the other and, our photo- 
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metric system has its own deficiencies. For example, because of tlie elongated stellar images 
on the CCD frames, it is hard to accurately define the center of the aperture, especially for 
the i?'-band images. This results in an extra error in the data reduction. Another drawback 
may be that the wavelength separation between the B' and R' bands is not large enough. 
Therefore, a photometric system with a light splitting device may be more suitable in search- 
ing for the optical lags. When achieving a temporal resolution of about or higher than 1 
minute, which is more than an order of rnagnitude shorter than the timescale of the op tical 



variability of S5 0716+714 fiSasada et al.ll2008l : iRani et all 12010a] : iChandra et al.ll201lD . the 



traditional quasi-simultaneous photometric system may also be able to detect the optical 
lags with high confidence level. 

We thank the anonymous referee for insightful comments and suggestions that helped 
to improve this paper very much. This work has been supported by the Chinese National 
Natural Science Foundation grants 10873016, and 11073032, and by the National Basic 
Research Program of China (973 Program), No. 2007CB815403. 
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Fig. 1. — An example frame taken with the novel photometric system. Labeled are target 
blazar and four reference stars used for flux calibration. The B\ V\ and R! "images" are 
marked for Star 1 for an example. North is to the up and east, to the left. 
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Fig. 2. — Light curves in B' , V, and R' bands for the entire monitoring period. For clarity, 
the B' and R' band hght curves are shifted by 0.2 and —0.5 magnitudes, respectively. 
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Fig. 3. — Intranight light curves on JD 2,454,088 in the B' , V, and R' bands. The associated 
small panels give the differential light curves of the check star (star 6) . 
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Fig. 4. — Intranight light curves on JD 2,454,090 in the B', V, and R' bands. The associated 
small panels give the differential light curves of the check star (star 6) . 
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Fig. 5. — Intranight light curves on JD 2,454,092 in the B', V, and R' bands. The associated 
small panels give the differential light curves of the check star (star 6) . 
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Fig. 6. — Intranight light curves on JD 2,454,097 in the B' , V, and R' bands. The associated 
small panels give the differential light curves of the check star (star 6) . 
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Fig. 7.— Color-magnitude diagrams for JDs 2,454,088, 2,454,090, 2,454,092, and 2,454,097. 
The left panel is for B' — V color, while the right is for B' — E! color. Different symbols 
represent data on different nights. They are fitted linearly and separately for each individual 
night and are shown with different lines. The solid lines are fittings to all points, with the 
corresponding correlation coefficients given at the upper-left corners. The errors are not 
plotted here for clarity. The mean magnitude error is 0.018 mags and the mean color error 
is 0.022 mags. 
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Fig. 8.— Color evolution during flares from 0.245 to 0.292 on JDs 2,454,088. The filled circle 
and square mark the beginning and ending points, respectively. The numbers indicate the 
points in time sequence. 
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Fig. 9. — ZDCF correlation and fitting results on JD 2,454,088. The dashed lines are 
Gaussian fittings to the points, and their peaks are marked with the vertical dotted lines. 
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Fig. 10. — ZDCF correlation and fitting results on JD 2,454,090. The dashed lines are 
Gaussian fittings to the points, and their peaks are marked with the vertical dotted lines. 
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Fig. 11. — ZDCF correlation and fitting results on JD 2,454,092. The dashed lines are 
Gaussian fittings to the points, and their peaks are marked with the vertical dotted lines. 
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Fig. 12. — ZDCF correlation and fitting results on JD 2,454,097. The dashed lines are 
Gaussian fittings to the points, and their peaks are marked with the vertical dotted lines. 



Table 1. Observational Log and Results in the B', V' and E! Bands 



Date 


Time 


Julian Date 


Exp. 


B' 


W 

err 


5B' 


v 


V 

' err 


5V' 


E! 


T3l 

-'''err 


5R! 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


2006 12 18 


16:03:37.0 


2454088.16918 


150 


14.340 


0.018 


-0.014 


13.745 


0.011 


0.008 


13.263 


0.017 


0.014 


2006 12 18 


16:06:25.0 


2454088.17112 


150 


14.298 


0.017 


-0.000 


13.740 


0.011 


-0.021 


13.271 


0.017 


-0.015 


2006 12 18 


16:09:14.0 


2454088.17308 


150 


14.325 


0.018 


-0.025 


13.751 


0.012 


-0.011 


13.293 


0.017 


-0.005 


2006 12 18 


16:12:02.0 


2454088.17502 


150 


14.320 


0.017 


-0.004 


13.738 


0.011 


-0.009 


13.292 


0.017 


-0.010 


2006 12 18 


16:14:51.0 


2454088.17698 


150 


14.316 


0.017 


0.011 


13.771 


0.011 


-0.011 


13.279 


0.017 


-0.001 



Note. — Table 1 is published in its entirety in the electronic edition of the Astronomical Journal. A portion is shown hbre 



for guidance regarding its form and content. i 
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Table 2. Variation Test Result 



JD 


1^1 




TP 


77) 

ryf 


77) 


771 

-ro.oi 


Var.'^ 


2,454,088 


26 


108 


23.705 


29.534 


25.166 


1.930 


Y 


2,454,090 


36 


148 


4.233 


16.408 


3.600 


1.763 


Y 


2,454,091 


37 


152 


1.599 


5.270 


2.817 


1.750 


Y? 


2,454,092 


36 


148 


8.342 


13.132 


8.296 


1.763 


Y 


2,454,093 


38 


156 


0.687 


2.584 


1.317 


1.740 


N 


2,454,094 


13 


56 


0.684 


0.523 


1.521 


2.465 


N 


2,454,097 


38 


156 


4.359 


11.649 


3.068 


1.740 


Y 



Table 3. Intranight Variation Amplitude 





Ab' 


Ay 




JD 


(mag) 


(mag) 


(mag) 


2,454,088 


0.262 


0.222 


0.233 


2,454,090 


0.208 


0.148 


0.152 


2,454,092 


0.182 


0.132 


0.172 


2,454,097 


0.206 


0.134 


0.110 



Table 4. Time Lags Between Different Wavelengths'^ 



JD 


B'- 


-V 


V- 


-R' 


B'- 


-R' 


ZDCF+GF 


FR/RSS 


ZDCF+GF 


FR/RSS 


ZDCF+GF 


FR/RSS 


2,454,088 


1.16±0.08 


1.26±1.92 


0.77±0.08 


0.09±2.04 


2.16±0.08 


2.40±2.16 


2,454,090 


5.48±0.51 


4.98±5.82 


25.35±0.81 


28.14±7.02 


26.46±0.92 


30.12±11.04 


2,454,092 


2.64±0.37 


2.52±3.48 


1.54±0.40 


3.30±2.94 


4.89±0.46 


6.78±4.26 


2,454,097 


-9.76±0.47 


-10.2±4.38 


4.60±0.40 


3.54±5.76 


-0.68±0.56 


-4.86±9.66 



The lags are in unit of minute. 



